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Abstract 

Interplanetary coronal mass ejections (ICMEs) can reach the Earth's 
magnetosphere causing magnetic disturbances. It can be measured by satel- 
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lite and ground-based magnetometers. Data from the ACE satellite and 
from the geomagnetic field was explored here via discrete wavelet transform 
(DWT). The increase of wavelet coefficient amplitudes of the solar wind pa- 
rameters and geomagnetic field data analysis are well-correlated with the 
arrival of the shock and sheath region. As an auxiliary tool to verify the 
disturbed magnetic fields identified by the DWT, we developed a new ap- 
proach called effectiveness wavelet coefficient (EWC) methodology. The first 
interpretation of the results suggests that DWT and EWC can be effectively 
used to characterize the fluctuations on the solar wind parameters and its 
contributions to the geomagnetic field. Further, these techniques could be 
implemented to real-time analysis for forecast space weather scenarios. 
Keywords: Wavelet analysis, Solar wind-magnetosphere interaction, 
Magnetogram data, Geomagnetic storm. 



i 1. Introduction 



2 Several years ago, in 1957, Hannes Alfven postulated that the solar wind 

3 was magnetized and that the solar-wind flow draped the magnetic field over 

4 comets, forming a long mag netic tail downstream in the antisolar direction 

5 ( iKivelton and Russell . Il995h . So, this solar plasma expands out from the 
e Sun driven by thermo-electrodynamical processes. The solar magnetic field 
7 propagates "frozen" in the solar wind in a spiral-like configuration due to the 
s Sun's rotation. The geomagnetic field can be considered as a dipole mag- 

9 netic field, limited by the changes in the solar wind density and velocity, 

10 and by the variation in the strength and orientation of the interplanetary 

11 magnetic field (IMF) ( iTasciond . Il988l ). The electrodynamical interaction be- 



2 



tween that solar plasma and the Earth's magnetized atmosphere generates a 
complicated interrelated current system such as magnetosphere current, tail 
current, ring current, field- aligned current and ionosp heric current, adding t o 



15 a lower electric circuit with atmospheric discharges (IMendes et al 



2005 
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In other words, the interaction between the solar wind and IMF with the 
magnetosphere-ionosphere-ground creates a great variety of co mplex pro- 



cesses, which generate geomagnetic activity (IMendes et all 120061 ) . 

During solar events, when the solar plasma in expansion incides upon 
the Earth's intrinsic magnetic field, a substantial transfer of energy into 
the terrestrial magnetosphere may take place. Then, the normally existing 
magnetospheric and ionospheric quiet currents are widened and intensified. 

Among other phenomena, one of the characteristic signature is the geo- 
magnetic storm, a depression in the horizontal component of the Earth's mag- 
netic field (H) at middle to low latitude. The key parameters that control the 
solar wind magnetospheric coupling are the direction, the strength and the 
duration of the interplanetary magnetic field (IMF). For example, intense 
magnetic storms (Dst < —100 nT) are caused b y IMF southward compo 



nent stronger than 10 nT at least for 3 hours (e.g. 



Gonzalez and Tsurutani 



19871 ). Solar wind speed and density also play an important role in the 



3i form ation of the ring current, though their exact role is still contrqversia 



32 (e.g., 



Huttunen et al. 



2005 



Gonzalez and Tsurutani 



Wang et al.l. 



1983). 



2003 



Fenrich and Luhmann . 



1998 



An important solar event is the coronal mass ejection (CME) because 
it can cause geomagnetic storms. They are observed near 1 AU and are 
called interplanetary coronal mass ejections (ICMEs). The term magnetic 
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37 cloud (MC) is used to characterize an ICME hav i ng a specific configuration 



of IMF and p 



Burlaga et al. 



asma density (e.g. 



Gosling . 



1990 



Klein and Burlaga , 



1982 



19811 ) . The MCs have values of plasma beta significantly lower 
than 1 which is the ratio of plasma pressure to magnetic pressure. Near 1 
AU, MCs have enormous radial sizes (0.28 AU), with an average duration 



43 solar wind speed 420 km / s (e.g ., 



2000 



Goldstein. 


1983; 


Lepping and Berdichevskv , 


In 


Goldstein 


(1983 


), it was suggested for the 



first time that MCs are force-free magnetic field configurations (that is, when 
VxB = a(r)B). The con stant a solution fo r a cylindrical symmetric force- 



free equation was given by 



Lundquist 



fll950h 



In this paper, we work with time series variations of solar wind parameters 
due to solar events to verify its behaviors before reaching the Earth, the 
"cause", and the "effect" of the solar wind-magnetosphere interaction, the 
geomagnetic storm, in the ground magnetograms. Several magnetic storms 
that occurred on April, 2001 have been analyzed. These magnetic storms 
were associated with magnet ic clouds occurrence a lthough shocks have also 

(Il993l ). magnetic clouds can be 



54 been observed. According to 



Farrugia et al. 



used for the study of the solar wind energy input to the magnetosphere since 
the IMF components vary smoothly with time and retain the polarity for 
relatively long intervals. The discrete wavelet transform (DWT) have been 
selected and used with three levels of decomposition in order to detect the 
small regularities in the solar parameters and the geomagnetic data. The 
rest of the paper is organized as follows: Section [2] is devoted to present 
a brief introduction of geomagnetic storms, showing an example of solar- 
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62 interplanetary-magnetosphere coupling. Section [31 the data and the analyzed 

63 period is presented. Section HI gives a discussion about the results. Finally, 

64 section [5] presents the conclusion of this work. 



65 2. Geomagnetics Storm 

ee The primary causes of geomagnetic storms are supposed to be strong 

67 dawn-to-dusk electric fields associated with the passage of southward directed 

6B interplanetary magnetic fields, Bs, passing the Earth for sufficiently long in- 

69 tervals of time. The solar wind energy transfer mechanism is magnetic re- 



conn ection between the IMF and the Earth's magnetic field ( IGonzalez et al 



71 11994 ). 

72 As consequence, the level of magnetosphere activity varies widely. Ge- 

73 omagnetic activity is classified by intensity and usually described by the 

74 variation of indices to distinguish between a quiet and an active day (oc- 

75 currence of storm or substorm). The index most used in order to quantify 

76 the effects on low latitudes is the Dst index, and recently, Sym-H. It repre- 

77 sents the variat ions of the H component due to changes of the ring current 



78 ( ITasciond . 11988) . The Sym-H is essentially the same as the traditional hourly 

79 Dst index. The main characteristic of the 1 minute time resolution Sym-H 
so index is that the solar wind dynamic pressure variation are more clearly seen 
si than indices with lower time resolution. Its calculation is based on magnetic 

82 data provided by eleven stations of low and medium latitude. Only 6 of 

83 the stations are used for its calculation of each month, some stations can be 

84 replaced by others depending on the data conditions. 

85 The principal defining property of a magnetic storm is the creation of 
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an enhanced ring current due to the increase of the trapped magnetospheric 
particle population. These particles present a drift due to magnetic field 
curvature and gradient, which leads the ions to move from midnight to dusk 
and electron s from midnight towar d dawn surrounding the Earth close to the 



90 dip equator (I Gonzalez et al. 



19941 ) 



By the time of the late recovery phase, when the ring current is symmet- 
ric and trapped on closed trajectories most of its energy has already been 
dissipated. Th e ring current becomes s ymmetric just after the minimum 



94 Dst is reached (IDaglis and Kozyra 



20021 ). The asymmetric structure (dur- 



ing the main phas e) and the symmetric str ucture (during the recovery phase) 



96 were observed by 



Ebihara and Ejiri 



( I2000l ) using data obtained by the polar 



orbiting satellite NOAA 12. 

It must be noticed that the quasi-trapped particles in the inner radiation 
belt can sink to the South Atlantic Magnetic Anomaly (SAMA) which is 



characte rized by a 



ioi tensity ( 



Pinto et al. 



Nishino et al. 



obal minimum in the Earth's total magnetic field in- 



19921 ). These particles can reach ionospheric heights. 



(120021 ) verified variation on the ionospheric parameters mea- 

103 sured by ionosonde at Cachoeira Paulista (Brazil). Also, the processes which 

104 cause energetic electrons to precipitate in the atmosphere are: magneto- 

105 spheric wave-particle interaction, lightning or artificially induced wave-particle 
loe interaction, drift-resonance in teractions and wave-partic le interactions gen- 
107 erated by plasma instabilities (jpinto and Gonzalez . 1989 ). In this paper, the 
los Vassouras station was chosen to study the geomagnetic variations over the 
109 Brazilian sector because taking into account somewhat the SAMA influence. 
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no 3. Dataset and Methodology 

in In this section, we will present the ACE satellites data used to characterize 

112 the variations on the solar parameters due to the ICMEs propagation. Also, 

in we will present the data used to analyze the effects of the magnetospheric 

n4 activity effects generated on the ground. Our study concerns to magnetic 

us events that occurred on April, 2001. These solar events occurred during the 

lie solar maximum of the 23th solar cycle. The method used here is based on 

ii7 the Discrete Wavelet Transform (DWT). It will be used to verify the increase 

us of wavelet coefficient amplitudes associated to the ICMEs propagation and 

n9 its effectiveness in a development of geomagnetic storms. 

120 3.1. ACE satellite data 

121 In this paper, the geomagnetic activity during the month of April, 2001 

122 was studied using ACE satellite dataset and Sy m-H index. Both of them are 



123 available at the NOAA web site (ISPIDR 



20081 ). 



124 Fig. Ufa) shows the behavior of the mean values of solar wind parame- 

125 ters and the Sym-H index. Each panel presents, from top to bottom, IMF 

126 components (Bx, By, Bz in GSE coordinate), plasma density, velocity (Vx 

127 component) and Sym-H index. Our interest is to characterize the variations 

128 on the solar parameters due to the ICMEs propagation. Also, we used the 

129 Sym-H index to verify the global magnetic field reductions during storms due 

130 to ICMEs propagation on April 2001. 

131 All time series represent ACE and geophysical parameters with 1 minute 

132 resolution. On April 2001, it happened several solar disturbances which were 

133 detected by ACE instruments. Some ICMEs that caused these geomagnetic 
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134 dis turbances were studied by IWu et al.l (120031 ) ; iHuttunen et al.l (120051 ) . 



Huttunen et al 



( 120051 ) used ACE satellites data to investigate which were 
the possible candidates to MCs events. In this paper, we will verify the in- 
crease of w avelet coefficient a mplit udes associated to the same ICMEs events 



studied by 



Huttunen et al 



139 (MCs) identified by 



(2005J1 Table [TJ shows three magnetic clouds 



Huttunen et al. 



(120051 ) in each one is presented the MC 
arrival and stop time, the inferred flux-rope type and Dst index (day and 
time (UT) of occurrence). The minimum Dst i ndex is important to char - 



U2 acterize the MC s geoeffectiveness. According to 



Zhang and Burlagal ( 119881 ); 



Bothmerl ( 120031 ) the geomagnetic response of a certain MC depends greatly 



on its flux- rope structure (Tabled! column 4). 



Table 1: Magnetic Clouds identified by ACE data on April, 2001 ([Huttunen et all 120051 ) . 
The columns from the left to the right give: Shock date (day and hour (UT)), MC start 
date (day and hour (UT)) ; MC end date (day and hour (UT)), inferred flux- rope type 
(e-g., 



Huttunen et al 



2005), the minimum value of the Dst index, it date of the previous 



one (day and hour (UT), if the sheath caused the storm, it is indicated by "sh") 



Shock 


MC, Start 


MC, Stop 


type 




Day and time (UT) of Dst m i n 


11, 15:18 


12, 10:00 


13, 06:00 


WNE 


sh(-259) 


11, 23:00 


21, 15:06 


21, 23:00 


22, 24:00 


WSE 


-103 


22, 15:00 


28, 04:31 


29, 00:00 


29, 13:00 


SEN 


-33 


29, 03:00 



145 3.2. Ground magnetics data 

we The magnetic stations considered in the analysis are Honolulu (HON), 

147 San Juan (SJG), Vassouras (VSS), Hermanus (HER) and Kakioka (KAK). 

us We are interested to examine how the effects within the magnetic storms on 

149 April, 2001 can affect the global space and time configuration of the ring 
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150 current at low-latitudes and the additional variations of the H-component 

151 related to the effect of equatorwards penetration of electric fields from the 



152 field-aligned current (IWu et al 



20041 ). These selected magnetic stations are 



153 used to calculated the Dst index, excluding VSS, because the Dst index 

154 reflects the magnetic variations related to changes of the ring current. We 

155 choose the ground magnetic observatory of VSS station because we are also 
we interested to understand better how these processes can effect the magnetic 
is? response to disturbed times in the SAMA. 

158 Another point of interest is the correlation of the response of these chosen 

159 stations due to the development of geomagnetic storms. The geographic and 
wo geomagnetic localization of these stations are shown in Table [2j 

Table 2: INTERMAGNET network of geomagnetic stations used in this study. 



Station 


IAGA code 


Geograph 


ic coord. 


Geomaj 


metic coord. 




Lat.(°) 


Long.(°) 


Lat.(°) 


Long.(°) 


Honolulu (United States) 


HON 


21.32 


-158.00 


21.59 


-89.70 


San Juan (Puerto Rico) 


SJG 


18.12 


-66.15 


27.93 


6.53 


Vassouras (Brazil) 


VSS 


-22.40 


-43.65 


-13.43 


27.06 


Hermanus (South Africa) 


HER 


-34.41 


19.23 


-33.89 


84.68 


Kakioka (Japan) 


KAK 


36.23 


140.18 


27.46 


-150.78 



Source: http:/ /wdc. kugi.kyoto-u.ac.jp/igrf/gggm/index. html (2010) 

161 Fig. |3]^a) shows the behavior of the H-component of the Earth's geomag- 

162 netic field and the Sym-H index detected on April 2001. Each panel presents, 

163 from top to bottom, the H-component obtained at HON, SJG, VSS, HER and 

164 KAK; and the Sym-H index. We could observe a magnetic signature of few 

165 geomagnetic storms with different magnitudes. The geomagnetic storms can 
lee consist of four p hases: sudden commen cement, initial phase, main phase and 



167 recovery phase (j Gonzalez et al.l . [1994J). The four major geomagnetic storm 
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168 which present all four phases are: a very-intense magnetic storm with mini- 

169 mum Dst = —271 nT at 24:00 UT on April 11, 2001, a second one started 

170 on April 18, 2001 with minimum Dst = — 114 nT at 08:00 UT, followed by 
m the third one started on April 21, 2001 with minimum Dst = —102 nT and 

172 a last storm occurred on April 28, 2001 with minimum Dst = —47 nT. 

173 3.3. Discrete Wavelet Transform (DWT) 

174 The discrete wavelet transform (DWT) has the following propriety: the 



larger amplitudes of the wavelet coefficients are associated with 
signal c hanges or "details" of hig her frequency. In the wor k of 



ocally abrupt 



Mendesetal 



(l2005al ) and the following work of iMendes da Costal (1201 ll ). a method for the 
detection of the transition region and the exactly location of this discontinu- 
ities due to geomagnetic storms was implemented. 

The DWT produces the so called wavelet coefficients at different levels 



i8i and it is proved that t heir amp 



182 of the analyzed data (IMallat 



itude s can be used to study the local regularity 



19971 ). As smaller is the amplitude as regular 



183 is the analyzed data. Therefore, where the amplitudes are large we can 



184 associate it to some disturbance on the signal (IMendes et al. 



2005ah 



The wavelet transform in level j + 1 is given by 



d{ +l = 2j29(m-2k)ci n 



(1) 



where g is a high-pass filter, di +1 is the wavelet coefficient at level J + 1, and 



c? m are the scale coefficients at level j. In this transform, 



4 +1 =2^Mm-2A;)c4 



(2) 
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185 and h is a low-pass filter. 

186 In this study, we considered j = as the most refine level of the multi- level 

187 decomposition which is associated to one minute data resolution. In other 
las word, c^, -0 is the mean time fluctuation computed from the raw dataset. 

189 In our case, the highest amplitudes of the wavelet coefficients indicate the 
wo singularities on the geomagnetic signal in association with the disturbed pe- 
rn riods. On the other hand, when the magnetosphere is under quiet conditions 

192 for the geomagnetic signal, the wavelet coefficients show very small ampli- 

193 tudes. In this work, we applied this methodology with Daubechies orthogonal 

194 wavelet function of order 2 (db2) with data sample rate of one minute time 

195 resolution. Also, we used with three levels of decomposition which are asso- 

196 dated to pseudo-periods of 3, 6 and 12 minutes. On the physical point of 

197 view, these perio ds are relat ed to the propagation of Alfven wave and the 

198 Pci 5 pulsations ( Sait J . 1969 ). 

199 In order to facilitated the visualization and the analysis of the wavelet 

200 coefficients, we developed a methodology called the effectiveness wavelet coef- 

201 ficients (EWC). This new method has an advantaged over the DWT method 

202 presented above. In the DWT method, we have to verified the wavelet de- 

203 composition on all the calculated chosen levels, in our case, three levels. Only 

204 the local regularities detected on the three levels of decomposition could be 

205 considered due to a physical process, in this ICMEs propagation. The 

206 EWC simplifies the DWT method by reducing all the wavelet decomposition 

207 levels for only one level. It corresponds to the weighted geometric mean of 

208 the square wavelet coefficients per hour. It is accomplished by weighting the 
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209 square wavelet coefficients means in each level of decomposition as following 

EWC = 4 SI, *! + 2 21, * + p) 

210 where N is equal to 60 because our time series has one minute resolution. 
2n The advantage of the EWC over the DWT is that it make easy the visu- 

212 alization and the analysis of the wavelet coefficients. In this new approach, 

213 we only have to analyze one decomposition level. The EWC is weighting 

214 according with the decomposition levels because the wavelet coefficients are 

215 calculated by filteri ng the origina l signal and by subsampling the resulting 



216 signal by a factor 2 (IMallatl . 119971 ). In this case, the first level of the decom- 

217 position has half of the number of the points of the original data, the second 

218 level of the decomposition has half of the number of the points of the first 

219 level of the decomposition and so on and so forth. 



220 4. Results and Discussion 

221 In this section, we first analyzed the solar wind parameters due to solar 

222 events that happened at April, 2001 to verify it behaviors before reaching the 

223 Earth, the "cause" of the geomagnetic storms. Second, we verified the "ef- 

224 feet" of the solar wind-magnetosphere interaction, the geomagnetic storm, in 

225 the ground magnetograms. Following, we examined the correlations between 

226 the "cause" and the "effect" . 

227 4.1. ACE 

228 Fig. [I] is divided in four panels. At top left, it shows the behavior of the 

229 solar wind parameters (IMF components (Bx, By, Bz in GSE coordinate), 
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230 plasma density and velocity (Vx)) and the Sym-H index for April, 2001. We 

231 highlighted the three MCs events shown in Table [1] using the dashed lines 

232 in red for shock time, in green for the MCs arrival time and in light blue 

233 for MCs stop time. At top right, at bottom left and at bottom right, the 

234 panels show the behavior of the square wavelet coefficients for these solar 

235 wind parameters at three first levels j = 1, 2 and 3, denoted by d 1 , d 2 , d 3 , 

236 respectively. The wavelet coefficients identified the variations of the ACE 

237 parameters due to solar wind events with different properties and degrees of 

238 disturbance. The MCs events shown on Table [1] were all identified by the 

239 increase of the wavelets coefficients at the three decomposition levels. 

240 The wavelet coefficients present an increase of amplitude during the first 

241 event shown in Table [TJ During the MC of WSE-type, the magnetic field 

242 vector rotates from the west (W) at the leading edge to the east (E) at the 

243 trailing edge, being south (S). Also, the Bz-component has the same sign 

244 during the MC propagation and it has the axis highly inclined to the eclip- 

245 tic (Unipolar MCs). The first MC event of WNE-type caused one extreme 

246 storm happened with a minimum Dst = -256 nT (shock) at 23:00 UT on 

247 11th of April. It was identified by wavelet coefficients at all three levels of 

248 decomposition during its shock time. A great number of wavelet coefficients 

249 present an increase of amplitude were found between the shock and the MC 

250 starting time at all three levels of decomposition. These wavelet structures 

251 are associated to the sheath, region of sola r wind bounded by the shock front 



252 and the ICME leading edge (jOwens et al. 



2005|). 
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Figure 1: ACE: (a) Magnetic field components (X, Y and Z), Plasma Density, Velocity (X-component) and Sym-H index; and 
(b,c,d) the wavelet coefficients {<P) 2 for j = f,2,3 and Sym-H index, respectively. 



The second MC event of SNE-type also presented an increase of wavelet 
coefficients amplitude. This event caused medium degree of disturbance of 
Dst index according to NOAA classification (Dst = - 33 nT). The plasma 
density presented wavelet structures during the shock and the sheath region 
at all three levels of decomposition. This could be a consequence of the low 
values of plasma density associated with this particular event in the solar 
wind. The third MC event with WSE flux rope type also was identified by 
wavelet coefficients between its shock and the MC starting time. Its wavelet 
structures were very similar to the first MC. 



Huttunen et al. 



(120051 ). the 



Apart of the three MCs events identified by 
wavelet coefficients present an increase of amplitude during other periods. 
This fact could be associated to the occurrence of ICMEs observed by the 

jarge Angle and Spectrometric Coro- 



Solar anc 
naerafic ( 


Heliospheric Observatorv/ 


Cane and Richarc 


son, 


2003) 


Cane and Richardson 


(2003 


) are listed 



estimated time of related disturbance in the upstream solar wind, the start, 
and the end times of the ICME, the quality of the estimated boundary times 
(1, most accurate; 3, ill-defined) and the minimum value of the geomag- 
netic Dst index respectively. In the last column, "2" indicates whether the 
ICME has been reported as a MC which can be modeled by a force-free flux 
rope < \http : / /lepmfi. gsfc.nasa.gov/mfi/mag-cloudjpubl.htmiy , "0" indi- 
cates that the field shows little evidence of rotation, i.e., the ICME is not 
reported as cloud. The other four ICMEs studied by them are not reported 
as MCs. However, the th ree MCs reported by them are the same of the work 
( 12005k 



of 



Huttunen et al. 
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278 Altogether, the wavelet structures was increased during the shock and 

279 the sheath region at all three levels of decomposition for all the MC events 

280 mentioned above. 

281 Fig. [2] shows EWCs on April, 2001 for the ACE dataset. On the verti- 

282 cal axis, each panel presents the the weighted geometric mean of the square 

283 wavelet coefficients per hour, in other words, the EWCs. And, on the hori- 

284 zontal axis, the time of day in Universal Time (UT). On the analyses of the 

285 IMF components, the wavelet coefficients presented higher amplitude during 

286 the shock and the sheath region of the first and third MC. However, the 

287 second MC was better detected by EWCs of plasma density than by the 

288 IMF components, during its shock and sheath region. Also, there was an in- 

289 crease of EWCs amplitude of all solar parameter during the third shock and 

290 sheath region. This characteristic depends of the interplanetary disturbance 

291 and the physical processes involved such as: magnetic reconnection, viscous 

292 dissipation (type interactions) or the mixture of both cases. 



2003). 



Table 3: ICMEs identified by LASCO data on April, 2001 (jCane and Richardson , 
The columns from the left to the right give: Shock date (day and hour (UT)), MC start 
date (day and hour (UT)), MC end date (day and hour (UT)), Quality, the minimum 
value of the Dst index (nT) and if the the cloud was identified. 



Shock 


ICME, Start 


ICME, Stop 


Quality 


Dst m i n 


MC? 


4, 14:55 


5, 11:00 


7, 03:00 


3 


-38 





8, 11:01 


8, 19:00 


10, 10:00 


2 


-51 





11, 13:43 


11, 22:00 


13, 07:00 


2 


-257 


2 


13, 07:34 


13, 09:00 


14, 12:00 


1 


-66 





18, 00:46 


18, 12:00 


20, 11:00 


2 


-100 





21, 16:01 


21, 23:00 


23, 08:00 


1 


-104 


2 


28, 04:31 


28, 14:00 


1, 02:00 


2 


-33 


2 



16 



-0 




Time (days) 



Figure 2: The effectiveness waveiet coefficients on April, 2001 for the ACE dataset. 



293 Also, it can be noticed an smaller EWCs on the IMF components that 

294 corresponded to the ICMEs propagation between the days of 4-7, 8-10, 11,13 

295 and 18-20, see Table [31 Once more, the larger EWC amplitudes of these solar 

296 parameter occur during the shock and sheath region. 

297 4.2. Ground magnetics data analysis 

298 Fig. [3] shows four panels presenting (a) the ground magnetograms and 

299 (b-d) the first three decomposition levels, respectively. Each panel displays, 

300 from top to bottom, the data or results for the magnetic stations of Honolulu 

301 (HON), San Juan (SJG), Vassouras (VSS), Hermanus (HER) and Kakioka 

302 (KAK), respectively, plus the Sym-H index. We observe that the increase on 

303 the coefficient amplitudes during the sudden storm commencement (SSC), 

304 and also during the main phase of geomagnetic storm. It is possible to no- 

305 tice that the wavelet coefficients presented an increase by the arrival time 

306 of the MCs when it reaches Earth's magnetosphere. The MC candidates 

307 are displayed on Table [T] and [3J All the magnetic stations showed similar 

308 pattern of magnetic behavior response to the storms with just few s i ngular- 
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fl2005ah 



309 ities in the three levels of wavelet decomposition. As 

310 suggested, this behavior may be related to differences in the magnetic coor- 

311 dinates, type of magnetometer used, local time, ground conductivity and Sq 

312 currents effects. 

313 In Fig. HI the first shock presented the higher EWCs, the same happened 

314 in Fig. [2j Similar to Fig. [2J the second shock was not well detected in Fig. HJ 

315 The third shock also presented the higher EWCs during the shock and its 

316 behavior was similar to Fig. [2j 
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Figure 3: Ground: (a) magnetograms and Sym-H index; and (b,c,d) the wavelet coefficients (d?) 2 for j = 1,2, 3 and Sym-H 
index, respectively. 
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Figure 4: The effectiveness wavelet coefficients on April, 2001 for the ground magnetograms. 
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Also, the EWCs shows an increase of values associated with the shocks 
listed on Table |3] and it was well-correlated to the arrival time of the ICMEs 
when it reaches Earth's magnetosphere. 



320 4-3- ACE and Ground data correlation analysis 

321 In the ACE parameter analysis, there was an increase of wavelet co- 

322 efficient amplitudes and EWCs of the IMF components during shock and 

323 sheath region. It can start a geomagnetic storm because during the pas- 

324 sage of a sheath, the shock compression, turbulence, magnetic field draping 

325 or shock he l iosph eric current sheet may lead to sout hward magnetic fields 



326 f lBadruddinl . 120021 ). As lGonzalez and Tsurutanil (119871 ) discussed, a geomag- 



327 netic storm occur when the interplanetary magnetic field turns southward 

328 and remains southward for an prolonged period of time. 

329 The ACE satellite is localized in the Lagrangian Point LI. On April, 2001, 

330 it could detected solar events about 30 to 80 minutes before they arrived to 

331 the Earth's magnetosphere (see Fig. If the DWT could be applied online 

332 then it could help to forecast a future geomagnetic disturbance. Generally, 

333 models that predict geomagnetic activity work well at minimum solar activ- 

334 ity Then, the capacity of the wavelet method to work correctly in high solar 

335 activity is an important advantage for its future use as a sophisticated space 

336 weather tool. 

337 In all the ground magnetograms analysis, the amplitude increase of wavelet 

338 coefficients were more frequent during the sudden storm commencement 

339 (SSC), and also during the main phase of geomagnetic storm. The SSC 

340 occurs due to sudden impulse due to the arrival of the interplanetary shock 

341 structure what generally coincides with the increase ram pressure (initial 

21 
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Figure 5: ACE data: Minutely time delay of the solar events arrival on the Earth's 
magnetosphere for April 2001. 

342 phase) followed by the decrease of geomagnetic field (main phase) which in- 

343 dicates sustained southward interplanetary fields in the sheath region or/and 



344 during the MC propagation (IGonzalez et al 



19941 ) . It was possible to notice 



345 a simultaneous increase of EWC values on Fig. |2]andHl respectively at, the 

346 shock and SSC regions; and the sheath and the main phase regions. Unfor- 

347 tunately, the weighting the square wavelet coefficients means by hour (EWC 

348 hourly values) disguises the arrival delay of the shock and sheath region at 

349 the Earth's magnetosphere. 

350 On the analyses of ground magnetograms, the highest coefficient am- 

351 plitude were coincident in time, showing that the whole magnetosphere is 

352 globally affect e d, at 
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east in the time resolution considered as discussed in 
(j2005al ). The amplitude of the EWCs are related to geoef- 

354 fectiveness of the solar events. It is also worth to mention that the small 

355 amplitudes of the wavelet coefficients mean that the energy transfer process 

356 is smooth; while the large amplitudes indicated that were an impulsive energy 
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injections superposed to the smooth background process. 

Table 14.31 shows the determination coefficient matrix between the ten 
dataset parameters (Bx, By, Bz, density, Vx, HON, SJG, VSS, HER and 
KAK) used in our study. The determination coefficient is defined as the 
square of the correlation coefficient. One of the reasons to use the coefficient 
of determination instead of the correlation is to compute the statistics in or- 
der to determine the size or magnitude of the relation between two variables. 
It is interpreted as the percentage of variability of the response variable ex- 
plained by the regression model. The correlation coefficient measures linear 
association. Though in spa ce geophysic s both the determination and corre- 



367 lation coefficient are used (IReiff 



19831 ). In our case, we prefer to use the 



368 determination coefficient due to its interpretation of linear regression. 



Table 4: Determination coefficient matrix between the ten parameters used in this study 
of solar wind and magnetogram datasets. 





B x 


By 


B z 


rip 


Vx 


HON 


SJG 


VSS 


HER 


KAK 




100. 




















By 


90.82 


100. 


















B z 


92.54 


90.82 


100. 
















n v 


18.66 


21.34 


25.40 


100. 














V x 


53.14 


33.29 


42.51 


0.57 


100. 












HON 


52.71 


58.37 


46.65 


10.50 


22.56 


100. 










SJG 


29.92 


29.37 


26.21 


4.84 


27.14 


83.91 


100. 








VSS 


17.81 


17.56 


15.29 


2.25 


21.34 


72.08 


97.02 


100. 






HER 


45.97 


48.02 


40.96 


8.41 


27.56 


95.45 


89.49 


79.92 


100. 




KAK 


40.45 


41.34 


35.64 


7.56 


28.94 


92.92 


95.84 


88.55 


97.22 


100. 



369 As we mentioned above, the larger values of determination coefficient 

370 between the five magnetic stations reinforces the statement that whole mag- 

371 netosphere is globally affected by the solar disturbances, see Table 14.31 At 
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372 the same time, VSS station presented the lower values of determination co- 

373 efficient between the five magnetic stations. This behavior of VSS can be 

374 explained by its localization which is in a region very peculiar, under the 

375 influence of the SAMA that is characterized by a glo bal minimum in the 

376 Earth 's total magnetic field intensity. As suggested by iPinto and Gonzalez 



377 ( 119891 ). the SAMA region can be compared to the auroral region mainly 

378 during geomagnetic storms due to the development of some current- driven 



379 plasma instab i 
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ities. This peculiar behavior of VSS was also observed by 



( I2005a[ ) in their study using several magnetic stations and ge- 



381 omagnetic storms. 

382 Also, it is possible to notice larger values of determination coefficient be- 

383 tween the Bx, By and Bz solar wind parameter is expected because their 

384 are the magnetic component of the same magnetic vector. However, between 

385 the solar wind parameters and the five magnetic stations the determination 

386 coefficients present lower values due to the amplitude of the wavelet coef- 

387 ficients are related to the propriety of the db2 wavelet to reproduce linear 

388 function locally, this means that it is sensible to local variations on the solar 

389 parameters. In other words, not every fluctuation of the solar parameters 

390 has enough energy to disturbed the Earth's magnetosphere but it is detected 

391 by DWT analysis. 



392 5. Conclusions 

393 In this work, we investigate the magnetic disturbances that occurred on 

394 April, 2001 sing a new methodology based in the amplitude of the discrete 

395 wavelet coefficients applied to solar wind and magnetogram datasets. We se- 
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396 lect the solar parameters detected by the ACE satellite and the H-component 

397 of magnetic field measured by ground magnetometers at HON, SJG, VSS, 

398 HER and KAK. The magnetic station of VSS were selected study to 

399 verify the geomagnetic responses to this event under the region of the SAMA 
wo influence. 

401 In response to the disturbances on April, 2001, the results show that the 

402 increase of wavelet coefficient amplitudes associated with the solar wind pa- 

403 rameters and H-component were well correlated with the shock and sheath 

404 regions. In the ground magnetograms analysis, the amplitude increase of 

405 wavelet coefficients were well correlated with the sudden storm commence- 

406 ment (SSC) and the main phase of geomagnetic storm. These phase of the 

407 geomagnetic are associated with the arrival of the shock and sheath regions 

408 on the Earth's magnetosphere, respectively. Our results also show that the 

409 VSS station presented a peculiar behavior that could be explained due to it 

410 is located under the influence of the SAMA. 

411 The effectiveness wavelet coefficient (EWC) methodology facilitates the 

412 visualization and the analysis of the wavelet coefficients because it only 

413 have one decomposition level. This new appro ach suggests a new r epre- 



sentation of the r e sults. The previous studies of 
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f l2005al ) and 

Mendes da Costal ( 1201ll ). employed three decomposition levels. Our method- 



416 ology could be used in a semi-automatic way to characterize the solar wind 

417 and ground magnetic disturbances. It could be important to improve the 

418 knowledge of the peculiarities and the effects produced by solar events on 

419 the Earth's magnetosphere and further it could be implemented to real-time 

420 analysis for forecast space weather scenarios. And also, it could enable us to 
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421 study certain aspects of processes involved in the solar wind-magnetosphere 

422 interaction. 
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